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ABSTRACT 

We present new results of a study of the sub-millimetre continuum emission from a 
sample of 9 radio galaxies and 4 radio-quiet quasars at redshifts z = 0.75 — 4.26. 
The observations were made at 800/im, using the single-element bolometer UKT14 
on the James Clerk Maxwell Telescope (JCMT), reaching a typical r.m.s. sensitiv- 
ity of <J rm s ~ 4 mjy and represent some of the deepest submillimetre extragalactic 
measurements made to date. Three detections were achieved, of which two are secure 
(4C41.17, Dunlop et al. (1994) and H1413+117) and one (53W002) is tentative, whilst 
comparable upper- limits were obtained for 7 of the 10 remaining sources. We use these 
data as the motivation for a detailed discussion of the conversion from submillimetre 
and millimetre continuum fluxes to dust /gas masses and star formation rates at high- 
redshift, and determine these quantities from our own and other data on high-redshift 
radio-galaxies and quasars. In particular we have investigated the impact of the four 
main sources of uncertainty in deriving physical quantities from such data, namely 
i) potential contamination by galactic cirrus, ii) uncertainty in the value of the dust 
rest- frequency mass-absorption coefficient, iii) difficulty in constraining the dust tem- 
perature, and iv) estimation of the appropriate gas:dust ratio in high-redshift objects. 
Our discussion emphasises how important it will be to quantify and, where possible, 
minimize such uncertainties (via, for example, appropriate observational strategies) 
in order to fully capitalize on the ten-fold improvement in sensitivity offered by the 
imminent arrival of the next generation of bolometer arrays, such as SCUBA on the 
JCMT. 

Taking these uncertainties into account we show that whilst the high-redshift 
galaxies discussed in this paper are amongst the most dust-rich and luminous objects 
discovered to date, their far-infrared properties are more comparable with those of 
the most luminous nearby interacting galaxies than with those expected of primaeval 
giant ellipticals. This conclusion is rather insensitive to the adopted dust temperature, 
and the appropriateness of our adopted gas:dust ratio is supported by the molecular 
line detections of lensed objects. Indeed, despite all of the uncertainties peculiar to 
studying galaxy evolution at sub-millimetre wavelengths, the current uncertainty over 
the true value of f2 is probably the most important factor affecting our conclusions. 
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1 INTRODUCTION 

Over the past four years, we have undertaken a number of 
sub-millimetre observations of high-redshift radio galaxies 
and quasars (using the single channel 3 He bolometer UKT14 
on the JCMT) aimed at determining the evolutionary sta- 
tus of these objects on the basis of the strength of their 
rest-frame far-infrared continuum emission. Since the arrival 
of the new generation of sub-millimetre bolometer arrays is 
now imminent, this seems an appropriate time to summa- 
rize the existing observational results (including our own 
data presented here), to review the main sources of uncer- 



tainty which afflict the meaningful physical interpretation of 
such data, and, in the light of recent Keck and HST results 
(Steidel et a£.1996, Illingworth 1997), to re-appraise the im- 
portance of deep sub-millimetre observations for the study 
of galaxy evolution and formation (Hughes 1996). 

The original motivation for attempting to detect sub- 
millimetre emission in high-redshift galaxies, despite the rel- 
ative insensitivity of existing sub-millimetre detectors was 
the realistic possibility that at least some massive galaxies 
formed the majority of their stars in a relatively short (< 1 
Gyr) starburst at high redshift, and that the resulting nec- 
essarily high star-formation rates (~ 100 — > WOOMQyr -1 ) 
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could give rise to strong rest-frame far-infrared thermal 
emission from dust. 

This motivation remains as strong as ever despite the 
recent discovery at optical wavelengths of a population of 
star-forming galaxies at redshifts z ~ 3 by Steidel et al. 
(1996). While the space-density of these Lyman- limit galax- 
ies at z ~ 3 indicates that we may be seeing the progenitors 
of a substantial fraction of present-day bright galaxies, their 
rather moderate star formation rates (~ ldM^yr^ 1 ) would 
have to be sustained for virtually a Hubble time to pro- 
duce the most massive elliptical galaxies which exist in the 
present-day Universe (of stellar mass ~ 10 12 Mq). The im- 
plication is that, while in the Lyman-limit galaxies we may 
be seeing the formation of the bulges of spiral galaxies and 
the cores of some ellipticals (an interpretation supported by 
their compact spheroidal appearance; Giavalisco et al. 1996), 
certainly for the most massive galaxies the star-formation 
rates must have been substantially greater (> lOOMoj/r -1 ) 
at higher redshifts. 

Indeed, we now possess rather strong evidence that 
this must have been the case. First, the basic properties 
of present-day elliptical galaxies, namely low molecular gas 
and dust masses, < 10 Mq (Lees et al. 1991, Knapp & Pat- 
ten 1991, Wilkind & Henkel 1995), enormous stellar masses 
~ 10 11 — 10 12 Mq, determined from the K-band luminosi- 
ties (S> 2 x L* , Taylor et al. 1996) of the host galaxies of 
low-z radio galaxies, RLQs and the most luminous RQQs, 
and uniform optical-IR colours (Bower, Lucey & Ellis 1992) 
that are dominated by a well-evolved stellar population, in- 
dicate that the bulk of their stars were formed in a relatively 
short-lived star-burst at high redshift. Second, rather more 
severe constraints on what the phrase 'high-redshift' must 
actually mean have been recently provided by Dunlop et 
al. (1996), who have presented spectroscopic evidence that 
star formation ceased in a z = 1.55 radio galaxy at least 
3.5 Gyr prior to the epoch at which the galaxy is observed. 
Such a large age formally excludes an Einstein-de-Sitter cos- 
mology for a Hubble constant Ho > 50 kms _1 Mpc _1 and 
thus indicates an extreme formation redshift (z > 5) for 
this galaxy in almost any cosmological model. This object is 
not a freak; a second, apparently even older radio galaxy has 
now been discovered at a similar redshift (Dunlop 1997) and 
spectroscopy of radio quiet ellipticals in the cluster around 
the 2 = 1.206 radio galaxy 3C324 indicate similarly old ages 
for the vast majority of their stellar populations (Dickinson 
1997). Third, a radio galaxy has recently been discovered at 
z = 4.41 (Rawlings et al. 1996); this demonstrates that at 
least some massive galaxies were in place at early epochs, 
and its properties (e.g. low rest-frame ultraviolet flux) are 
not obviously those of a young star-forming system. In sum- 
mary, despite the important discovery of the star-forming 
population at z ~ 3 (Steidel et oZ.1996), a growing body 
of evidence in fact indicates that we have yet to discover 
the formation epoch of the most massive galaxies and that, 
given the available cosmological time, the formation of these 
objects must have been shortlived and hence should be spec- 
tacular in at least one wavelength regime. 

These last two studies illustrate the key role which radio 
galaxies continue to play in studies of high-redshift galaxy 
evolution and formation. The main reason for this is that all 
low-redshift radio galaxies are associated with giant ellipti- 
cal galaxies (Owen & Laing 1989), and thus it is reasonable 



to assume that they can be used to trace the evolution of 
the most massive galaxies back to high redshifts, and early 
cosmic epochs. 

However, while we can be reasonably confident that the 
hosts of high-redshift radio sources are the progenitors of gi- 
ant elliptical galaxies, recent years have seen a growing con- 
cern that, at least in the most powerful sources, the direct 
or indirect effects of their active nuclei might stymie any 
attempt to determine their evolutionary status at optical- 
infrared wavelengths. For example, the discovery of multi- 
modal optical structures aligned along the radio axis of high- 
redshift radio galaxies (Chambers et al. 1987; McCarthy et 
al. 1987) seemed at first sight to indicate the assembly of gi- 
ant elliptical galaxies from the merging of lower-mass clumps 
at relatively recent redshifts, as expected in certain hier- 
archical models of structure formation. However, in subse- 
quent years there has been much debate about whether this 
is true, and about the evolutionary status of high-redshift 
radio galaxies in general. Early arguments that z > 2 radio 
galaxies contained old evolved stars (e.g. Lilly 1988; Lilly 
1989) have been questioned (e.g Chambers & Chariot 1990; 
Eales & Rawlings 1993) but there has been a distinct lack 
of proof that any of the stellar systems are demonstrably 
young (although arguments have been made that this is 
the case; e.g. 53W002 - Windhorst et al. 1991; 0902+34 
- Eales et al. 1993). The few hints of 'primaevality', gleaned 
from ultraviolet-infrared observations, have remained hard 
to prove largely because of the way in which high-redshift 
radio galaxies are typically selected by virtue of their ex- 
treme radio luminosity. An inescapable consequence of high 
radio luminosity appears to be significant non-stellar emis- 
sion across the ultraviolet, optical and infrared bands (e.g. 
Rawlings & Saunders 1991; Dunlop & Peacock 1993; Eales 
& Rawlings 1996). In other words huge narrow emission 
lines, large ultraviolet luminosities and multi-modal struc- 
tures (which naively one might interpret as evidence for ex- 
tremely high star-formation rates and dynamically young 
systems) are now often attributed to to the action of the 
powerful jets or quasar light emanating from the active nu- 
clei of these extreme objects. On the other hand, even when 
radio galaxies are faint (e.g. the z = 4.41 radio galaxy stud- 
ied by Rawlings et al. 1996) it is impossible to disprove huge 
star-formation rates since, as is the case for some local ultra- 
luminous IRAS galaxies (Sanders et al. 1988), star formation 
activity may be enshrouded in dust. 

This latter point provides the prime motivation for 
the study described here. A more robust signature of mas- 
sive star formation than ultraviolet continuum and emission 
lines is intense FIR emission from dusty, molecular mate- 
rial, where the rate of dust production is proportional to 
the star formation rate. The dust is heated primarily by the 
embedded O and B stars which evolve quickly and disperse 
their surrounding material on similarly short timescales 
(~ 10 7 yrs., Wang 1991). Hence the FIR luminosity provides 
a measure of the current formation rate (SFR) of massive 
stars, 

SFR = *10- 10 ^f^ M^yr- 1 (1) 

where * = 0.8-2.1 (Scoville & Young 1983, Thronson 
& Telesco 1986). If galaxies at high redshift have FIR lu- 
minosities comparable to or greater than low-z ULIRGS 
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Figure 1. The anticipated evolution of the UV to millimetre 
wavelength rest-frame spectral energy distribution of elliptical 
galaxies. Models corresponding to galactic ages between 2-15 
Gyr are shown, assuming a Salpeter IMF with a lower mass limit, 
mi = O.OIMq. The SFR (il>(t) = (m gas /m ga i)ipo) is proportional 
to the fractional mass of gas in the galaxy and assumes an in- 
tial SFR Vo = 100M© /yr. The predicted SEDs are taken from 
Mazzei, de Zotti & Xu (1994). 

(L FIR /L Q > 10 12 ), hence SFRs > lOOMoyr" 1 , then it is 
possible that they convert 10 11 — 1O 12 M0 of gas into stars 
in a burst of duration < 1 Gyr. 

In the Milky Way and local disc galaxies a significant 
fraction (30%) of the bolometric luminosity (L bo i) is re- 
radiated at FIR wavelengths, and hence the SFR (Miller & 
Scalo 1979, Kennicutt 1983) and the ratio Lfir/ 'Lbol cannot 
have evolved much with look-back time. However the situa- 
tion is very different for elliptical galaxies. Mazzei, de Zotti 
& Xu (1994) have modelled the photometric evolution of el- 
liptical galaxies and show that, whilst at the current epoch 
ellipticals emit < 1% of their bolometric luminosity at FIR 
wavelengths, within the first 1 — 2 Gyr of their formation this 
fraction was significantly higher with Lfir/ L bol ~ 0.3 (see 
figure 1). 

Whilst the details of the evolution are sensitive to the 
assumed initial mass function (IMF) and SFR, (where the 
steeper IMF and higher SFR produces a more luminous, 
but shorter, burst of starformation) , it is hard to escape the 
conclusion that the formation of giant elliptical galaxies is 
expected to be a spectacular and luminous phenomenon in 
the rest-frame far-infrared. Indeed this should be true ir- 
respective of whether ellipticals form via the collapse of a 
single gaseous halo, or grow through the rapid merging of 
smaller mass clumps at high-z, particularly since the mass 
dependence of metallicity and massdight ratio in elliptical 
galaxies indicates that early star-formation in massive el- 
lipticals should be strongly biassed towards high-mass stars 
(Zepf & Silk 1996). 

Any attempt to detect this thermal radiation from dust 
in high-redshift radio galaxies would fail were it not for 
the fact that the expected far-infrared emission peak (at 
A ~ 60 — 100/itm), due to grains radiating at tempera- 



tures of 30-70 K, is shifted into the submillimetre region 
(A > 300/im) at the redshifts of interest. Figure 2a demon- 
strates that, in an Einstein-de Sitter universe and at wave- 
lengths > 800/im, the effects of cosmological dimming on 
a typical starburst galaxy spectrum are offset between red- 
shifts z = 1 — > 10 by the strongly negative k-correction 
(which arises as the submillimetre filters effectively climb 
the Rayleigh- Jeans tail of the thermal dust emission with 
increasing redshift). The situation for a low density universe 
is not as advantageous (figure 2b) , where the observed flux 
density for a redshifted object continues to fall, albeit more 
slowly, in all submillimetre and millimetre wavelength pass- 
bands. 

In this paper we report on an attempt to exploit the 'de- 
tectability' of thermal emission from dust at high redshift 
to determine the evolutionary status of known high-redshift 
radio galaxies and quasars. The aim is to use sub-millimetre 
continuum photometry to estimate or at least constrain the 
mass of dust in high-redshift objects, with a view to deter- 
mining both the star formation rate and, more importantly, 
the mass of gas which has yet to be turned into stars in these 
potentially young galaxies. 

While gas masses can in principle be derived from 
molecular line observations, we have chosen to concentrate 
on a sub-millimetre continuum approach for three reasons. 
First, comparable assumptions and uncertainties also com- 
plicate attempts to derive an accurate H2 mass from de- 
tections of CO emission. Second, despite determined efforts, 
with the exception of BR1202-0725 (Omont et al. 1996, 
Ohta et al. 1996), no significant detection of molecular gas 
(through, e.g. CO line transitions) has been achieved in any 
high-z radio-galaxy or quasar (Evans et al. 1996, van Ojik 
et al. 1997, Barvainis & Antonucci 1996) unless it has been 
greatly magnified by gravitational lensing (Barvainis et al. 
1994, Eisenhardt et al. 1996). Third it is in sub- millimetre 
continuum astronomy that the greatest advances in sensi- 
tivity are likely to be forthcoming in the immediate future. 

As explained above, we have chosen to concentrate on 
radio galaxies because of their likely association with the 
high-redshift counterparts of massive elliptical galaxies, but 
we have also observed a few radio-quiet quasars in an at- 
tempt to confirm previously published detections. Our sub- 
millimetre observations are described in Section 2. In section 
3, we present a detailed point-by-point analysis of the var- 
ious uncertainties which afflict attempts to accurately es- 
timate the dust mass of a high-redshift galaxy from sub- 
millimetre observations. Then in section 4 we derive our best 
estimates of the dust masses of these high-redshift objects 
and use this information to assess their evolutionary sta- 
tus, by estimation of both the star formation rate, and the 
amount of molecular gas which, at the time of observation, 
has yet to be converted into stars. Finally, in section 5 we 
discuss explicitly, and endeavour to quantify the realistic un- 
certainties in these extrapolations, in particular the gas:dust 
ratio at high-redshift. 

Throughout the paper we assume qo = 0.5 and Ho = 
50kms _1 Mpc -1 unless stated otherwise, and correct pre- 
viously published data to the same cosmology to allow a 
comparison of all physical quantities. 
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Figure 2. Dependence of the observed flux density on redshift of the starburst galaxy M82 (Hughes et al. 1994) a t 350fim (dashed line), 
800/im. (solid line), 1100/im. (dashed-dotted line) and 2000^m (dotted line) for Ho = 1 (left-hand plot) and Qq = (right-hand plot) 
assuming Hg = 50 km s _1 Mpc -1 . 



2 SUBMILLIMETRE OBSERVATIONS 
2.1 Target selection strategy 

As figure 2 illustrates, if one wishes to attempt to detect 
high- redshift (z > 1) objects at sub-millimetre wavelengths, 
little if anything is lost in terms of sensitivity between z ~ 1 
and z ~ 10. Consequently, since the star-formation rate in 
elliptical galaxies is expected to be highest at early times 
(figure 1) it makes sense to start with the most distant 
known objects. We thus selected our primary targets (see 
Table 1) from among the most distant known radio galaxies, 
concentrating on objects which have been well-studied in 
the optical-infrared while also attempting to span a mean- 
ingful range of radio luminosities and redshifts. 4C41.17 and 
8C1435+643 are extremely luminous radio sources and each, 
at the time they were observed, held the title of the most 
distant known galaxy in the Universe. B20902+34 is an ex- 
tremely well-studied object and has been hailed as one of the 
best known candidates for a primeval galaxy (Eales et al. 
1993). 6C0032+412 and MG2141+192 arc two of the most 
distant known galaxies at more moderate radio luminosities, 
while 53W002 is the most distant mjy radio galaxy (Wind- 
horst et al. 1991). 3C257 was selected because it is the most 
distant 3CR radio galaxy and because near-infrared spec- 
troscopy has indicated that Lyman-a emission in this source 
is attenuated by dust. Two additional radio galaxies at in- 
termediate redshifts (z ~ 1) were selected because of cir- 
cumstantial evidence suggesting the presence of dust; 3C318 
is the most distant radio galaxy detected by IRAS (Heck- 



man et ai.1994) and 3C65 has an extremely red r-K colour 
(Lilly 1989, Dunlop & Peacock 1993). Lastly we observed 4 
RQQs, again spanning the redshift range (z = 1 — » 4), both 
to attempt to confirm or refute claimed IRAM detections at 
millimetre wavelengths, and to enable us to at least make a 
preliminary attempt to compare the far-infrared properties 
of radio- loud and radio-quiet objects at high redshift. 

2.2 Selection of observing wavelength 

Obviously we have no a priori knowledge of the exact spec- 
tral shape of the continuum at FIR-mm wavelengths in the 
rest frame of the high-redshift radio galaxies. However, since 
our aim is to test whether these galaxies are undergoing an 
intense burst of starformation at early epochs, it seems rea- 
sonable to assume that their rest-frame FIR-mm continuum 
spectra will be similar in form to those of massive galac- 
tic starforming regions and low-z starburst galaxies. At the 
highest redshifts the FIR spectral peaks of such starforming 
regions fall conveniently into the region of two submillime- 
tre wavelength atmospheric windows at 350/im and 450/im 
and it might seem that observations would be most prof- 
itably made at these wavelengths. However the 750/xm and 
850/im atmospheric windows provide more suitable and, in 
practice, preferential alternatives. At redshifts z > 2 the 
observed flux density ratio S^oo^m / Ssoo^m of a starburst 
galaxy spectrum is significantly smaller than that at z ~ 
(Figure 3) and is no longer sufficient to offset the observa- 
tional disadvantages of increased sky noise and reduced sky 
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Table 1. 800 fim photometry of high-redshift radio galaxies (RGs) and radio-quiet quasars (RQQs). The flux density limits quoted in 
column 6 are given at the 3<r level. The references in column 7 give the positions, rcdshifts and details of other relevant observations. 



Source name type z R.A. (J2000) Dec (J2000) SsoO/am (mjy) refs 



3C 318 


RG 





752 


15 


20 


05 


49 


+20 


16 


05.1 


<33 


1 


3C 65 


RG 


1 


176 


02 


23 


43 


48 


+40 


00 


52.7 


<11 


1,2 


PG1634+706 


RQQ 


1 


334 


16 


34 


28 


97 


+70 


31 


32.4 


<47 


3 


53W002 


RG 


2 


390 


17 


14 


14 


78 


+50 


15 


30.4 


6.9 ± 2.3 


4,5 


3C 257 


RG 


2 


474 


11 


23 


09 


40 


+05 


30 


17.8 


<11 


6,7,8 


H1413+117 


RQQ 


2 


546 


14 


15 


46 


26 


+11 


29 


43.7 


66 + 7 


9 


2132+0126 


RQQ 


3 


194 


21 


35 


10 


61 


+01 


39 


31.3 


<12 


10 


B2 0902+34 


RG 


3 


391 


09 


05 


30 


10 


+34 


07 


57.3 


<14 


6,11,12 


MG 2141+192 


RG 


3 


594 


21 


44 


07 


52 


+ 19 


29 


14.2 


<11 


6,8 


0345+0130 


RQQ 


3 


638 


03 


48 


02 


29 


+01 


39 


18.4 


<25 


10 


6C 0032+412 


RG 


3 


665 


00 


34 


53 


09 


+41 


31 


31.5 


<14 


8,13 


4C 41.17 


RG 


3 


800 


06 


50 


52 


36 


+41 


30 


31.5 


17.4+ 3.1 


14,15,6 


8C1435+643 


RG 


4 


252 


14 


36 


37 


19 


+63 


19 


14.2 


<13 


16,17 



1. Laing et al. (1983) 2. Stockton et al. (1995) 3. Schmidt & Green (1983) 4. Windhorst et al. (1991), 5. Windhorst et al. (1994), 6. 
Eales & Rawlings (1993), 7. McCarthy et a/.(1995), 8. Eales & Rawlings (1996), 9. Magain et al. (1988) 10. Schneider et al. (1991) 11. 
Lilly (1988), 12. Eiscnhardt & Dickinson (1992), 13. Rawlings et al. in prep, 14. Chambers et al. (1990), 15. Miley et al (1992). 16. 

Lacy et al. (1994) 17. Spinrad et al. (1995) 
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Figure 3. The predicted redshift dependence of the observed 
flux-density ratio S 400 nm I 'S 'goo |im calculated from the isothermal 
dust model fitted to the FIR-mm continuum data of M82 (Hughes 
et al. 1994). 



transparency, A trans , at 450/im compared to 800/im, under 
even the dryest atmospheric conditions (p.w.v. ~ 0.5 mm) 
when A irans (450pm/800>m) is < 0.7. 

Therefore all observations were initally made in the 
850pim and 750/xm windows using a broad-band filter at 
~ 800/xm (A c = 781^m, AA ~ 200/im). In the event of 
a detection at 800/j.m follow-up observations were made at 
450^im. An additional theoretical point also strongly favours 
the use of a filter at ~ 800/xm over shorter wavelength filters 
since, as discussed in §3.3.2, a single photometric measure- 



ment at A <C 800/im is virtually useless for constraining the 
dust masses in sources at z > 1. 



2.3 Observations and Results 

Sensitive continuum observations of radio galaxies and 
RQQs in the redshift range 0.8 < 2 < 4.3 were made 
at 800/im using the single-element 3 He bolometer UKT14 
(Duncan et al. 1990) on the JCMT during a series of 
runs between April 1993 and February 1996. Data was 
taken only under excellent submillimetre observing condi- 
tions (r800Mm — 0.35). A 65mm focal-plane aperture was 
used for all observations and the resulting beamsize of 16.5 
arcsec corresponds to physical diameters of ~ 140 — > 100 kpc 
at z = 1 — * 4. An azimuthal chop throw of 60 arcsec at a 
frequency of 7.81 Hz was used to subtract the sky emission. 
For each source a number of repeated, short (10-20 minute) 
observations were concatenated using a method described 
by Hughes et al. (1993). Total on-source integration times 
were typically 3-5 hours. Primary calibration was performed 
against Uranus and Mars, with secondary calibration against 
a variety of AGB stars and compact HII regions (Sandell 
1994). The data collected to date with UKT14 are sum- 
marized in Table 1. Note that, after the successful commis- 
sioning of the 0.1K bolometer array SCUBA (> 1997), one 
can expect typical 3a sensitivities of 12mJy and 1 mjy at 
450/im and 800/ira respectively with an on-source integra- 
tion of 3 hours. 



2.3.1 High-redshift radio galaxies 

A highlight of this study has been the clear detection of 
4C41.17 at 800/im, together with a significant upper limit 
of 3cr < 56mJy at 450/im (Dunlop et al. 1994). This de- 
tection was subsequently confirmed by an IRAM observa- 
tion at 1.25mm (Chini & Kriigel 1994) at a level consistent 
with the thermal dust model presented by Dunlop et al. 
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(1994). We have also achieved a marginal 3a detection at 
800/!?r!, of 53W002 at a considerably fainter flux density of 
6.9 ± 2.3 mjy. Yamada et al. (1995) have a similarly ten- 
tative detection of 12 CO (1-0) in 53W002 leading to an H 2 
gas mass of 2 x 10 12 Mq. Whilst we have had no opportu- 
nity to confirm our continuum result, treating the 800/im 
observation as a limit of 3a < 7 mjy (M d = 1.5 x 10 8 M Q , 
see §3.2) allows us to place an unlikely lower limit on the 
gas-to-dust ratio of > 13000 and hence casts doubt on the 
validity of the CO detection. At the time of observation 
it seemed somewhat surprising that no 800/oim detection of 
B20902+34 was obtained given the claimed IRAM detection 
at 1.3 mm (Chini & Kriigel 1996). However recent 105 GHz 
observations (Yun & Scoville 1996, Downes et al. 1996) now 
suggest that a significant fraction (probably the majority) of 
the 1.3 mm continuum in B20902+34 is due to non-thermal 
emission (see §3.1). No detections were achieved for the re- 
maining six radio galaxies at z> 2; our non-detection of 
8C1435+643 could be regarded as casting some doubt on 
the validity of the IRAM detection by Ivison (1995), but 
in fact our 3 sigma upper limit for this source is still just 
consistent with a grey-body fitted through the 1.25mm data 
point (see section 3.1). 

2.3.2 High redshift radio-quiet quasars. 

Continuum observations at 800/im (Barvainis et al. 1992; 
Isaak et al. 1994) and at 1.25 mm (Andreani, La Franca & 
Cristiani 1993, McMahon et al. 1994, Ivison 1995, Omont 
et al. 1996) have suggested that thermal emission from dust 
has been detected in high-redshift radio-quiet quasars. In 
an attempt to confirm some of these more marginal re- 
sults and facilitate direct comparison with our sample of 
high-redshift radio galaxies, we made 800/itm observations 
of two (2132+0126 and 0345+0130) of the three radio-quiet 
quasars detected by Andreani et al. (1993). In 2132+0126 
the 1.25 mm detection (11.7 mjy) and our measured 800/im 
flux density limit (< 12mJy) produce a spectral index 
a < 0.06 (£„ oc v a ) over the rest-frame wavelength range 
298/oim — ► 191/im. This is inconsistent with that expected 
from a thermal dust spectrum with a grain temperature 
T > 20 K. Therefore our result casts serious doubt either 
on the 1.25 mm IRAM detections or on the proposed ther- 
mal nature of the sub- millimetre continuum, unless the dust 
grains are significantly colder than observed in low-redshift 
quasars (Chini et al. 1989, Hughes et ai.1993). Our flux limit 
of 3a < 25mJy at 800/im for 0345+0130 is just consistent 
with the flux density expected from an extrapolation of the 
1.25 mm detection (assuming an isothermal 50 K spectrum 
with a grain emissivity index /3 = 2) and thus provides no 
additional constraint on the nature of the FIR spectrum or 
the dust mass in this object. 

In the light of the recent detections of the Cloverleaf 
quasar (H1413+117) at 1.25mm, 100/im and 60/^m (Bar- 
vainis et al. 1995), which indicated it to have a peculiar 
SED (figure 4), we dearchived and recalibrated the original 
submillimetre data of Barvainis et ai.( 1992). Our recalibra- 
tion at 350/xm and 450/im agrees with those of Barvainis 
et al. (1992), whilst suggesting a flux density at 800/im of 
66 ± 12 mjy, 50% higher than that in the original paper. 
This prompted us to make a new 800/tm observation of 
H1413+117 in March 1996 which confirmed our recalibra- 
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Figure 4. The mid-IR to millimetre spectral energy distribution 
of the Cloverleaf quasar (H1413+117). The open circles show the 
data presented by Barvainis et al. (1995). The solid circle shows 
our new observation at 800/jm that is consistent with various 
optically-thin isothermal greybody spectra (30 K, /3 = 2 - solid 
curve; 50 K, j3 = 1 - dot-dash; 76 K, = 1 - dot-dot-dot-dash). 
Also shown for completeness are greybodies assuming 30 K, fi = 1 
(dashed curve) and 50 K, (3 = 2 - dotted curve. The alternative 
greybody spectra are normalised at 1300/tm. 

tion with a detection of 66 ± 9mJy (R. Ivison & J.Stevens, 
priv comm.), giving a weighted mean of 66 ± 7mJy for the 
two independent data sets which differs from the published 
800/im detection of Barvainis et al. (1992) at the 2a level 
and is much more consistent with that expected from a grey- 
body fit to the other data-points (figure 4). The SED of the 
Cloverleaf is now one of the most well-defined for a high-z ob- 
ject in the FIR-millimetre wavelength regime and illustrates 
the difficulty in constraining the dust temperature within 
the range 30 K < T d < 70 K (figure.4, see also §3.3.2). 

2.4 Comparison of high-z and low-z AGN at 
submillimetre wavelengths 

In Table 2 we list all galaxies and radio-quiet quasars with 
redshifts z > 2 that have published detections at one or 
more wavelengths between 350/im — 1300/im, regardless of 
whether they have been confirmed. For the purposes of this 
paper a 'detection' is defined as any published photometry 
with a signal-to-noise ratio > 3. The reliability of any de- 
tection will be discussed when appropriate. 

In figure 5 we compare the 800/im flux densities and 
limits for the high-redshift galaxies listed in Table 2 with 
the 800/oim flux densities observed for a variety of starburst 
galaxies and AGN at low-redshift (Hughes, Davies & Ward 
1997). The redshift dependence of the 800/im flux density 
of the starburst galaxy M82, previously shown in fig. 2, is 
reproduced here, together with the dependence for multi- 
ples of xlO, xlOO and X1000 the rest-frame FIR luminosity 
of M82 (3 x 10 10 L Q ). This plot illustrates a number of ba- 
sic but important points in a rather transparent manner. 
First, with the exception of IRAS10214+4624 for which, in 
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this plot, we indicate the inferred 800pm flux density after 
correcting for lensing (Eisenhardt et al. 1996), there is an 
obvious absence of sources, at any redshift, with Ssoojxm < 7 
mjy. This simply reflects the current sensitivity limits of 
sub-millimetre bolometers such as UKT14; it is clear from 
the figure that our 800/im observations of the high-z sources 
are amongst the faintest 800/u.m observations made to date. 
Second, the lines on this plot indicate that for sources at 
z > 1, such flux density limits correspond to FIR lumi- 
nosities which are several hundred times greater than low-z 
AGN. Third, this plot shows that with an improvement in 
sensitivity of a factor of ten (as is promised by new submil- 
limetre bolometer array detectors such as SCUBA on the 
JCMT) it should be possible to detect objects comparable 
to the most luminous starburst galaxies and ULIRGS seen 
in the local universe (e.g. Arp220 and Mrk231), with lumi- 
nosities Lfir > 10 12 Lq out to z ~ 10. We also note that 
little effort has been invested in making submillimetre ob- 
servations of galaxies in the intermediate redshift interval 
z = 0.1 -> 1. 

The important physical property which can be derived 
directly from these sub-millimetre detections/upper-limits 
is the mass of dust in a given galaxy. From this one can 
then go on to assess the evolutionary status of the object 
in question, either through estimation of the 'current' star- 
formation rate, or by estimating the gas mass of the object 
(by adopting a given gas:dust ratio). However, in several 
recent papers this crucial process has been described with 
misleading simplicity. In fact, even the first stage in this pro- 
cess (calculating the dust mass) has to be undertaken with 
great care because there exist at least five potential sources 
of significant uncertainty. Thus, before attempting to derive 
the physical properties of the objects listed in Table 2, to 
assess their evolutionary status (which we do in section 4), 
in the next section we present a detailed description of these 
5 sources of uncertainty, assess their relative importance and 
outline the prospects for improvements in each area. 



3 UNCERTAINTIES IN DETERMINING THE 
DUST MASSES OF HIGH-REDSHIFT 
OBJECTS 

In this section we describe, roughly in order of increasing 
subtlety, the problems which afflict the accurate determi- 
nation of dust masses from sub-millimetre observations of 
high-redshift galaxies. First we consider how one can estab- 
lish that the detected emission is indeed from dust, rather 
than being, for example, synchrotron emission (an impor- 
tant issue for our radio galaxy targets). Next we address the 
issue of how confident one can be that it is the high-redshift 
target that has been detected, and not merely foreground 
galactic cirrus. Then, assuming that it can be established 
that dust emission from the high-redshift target has indeed 
been detected we consider the effect of uncertainties both 
in the intrinsic properties of the dust and in its tempera- 
ture, and discuss what can be done to minimize these un- 
certainties. Lastly, we briefly highlight the often overlooked 
implications of the current uncertainty in the values of cos- 
mological parameters. 



3.1 Potential Problem 1: Have we detected 
emission from dust? 

To prove that any sub-millimetre emission detected from 
high-redshift galaxies is due to dust one must ideally demon- 
strate that the sub-millimetre SED is rising too steeply to 
be due to self-absorbed synchrotron radiation (e.g. Chini et 
al. 1989a, Hughes et al. 1993). In practice this is extremely 
difficult to achieve (both due to achievable signal-to-noise, 
and the fact that at high-redshift one is not always observing 
the Rayleigh- Jeans tail of the expected thermal emission), 
but at the very least it is necessary to demonstrate that the 
observed sub-millimetre spectrum is rising towards shorter 
wavelengths. 

A direct determination of a positive sub-millimetre 
index a (where f v oc v a ) at high redshift is lim- 
ited to submillimetre/millimetre continuum observations 
of only 4 objects; H1413+117 (Barvainis et al. 1992, 
1995), IRASF10214+4724 (Rowan-Robinson et al. 1993), 
4C41.17 (Dunlop et al. 1994, Chini & Kriigel 1994) and 
BR1202-0725 (Isaak et al. 1994, McMahon et al. 1994). 
However, in certain circumstances, even on the basis of a sin- 
gle sub-millimetre or millimetre wavelength, one can be rea- 
sonably confident that the detected emission lies above any 
reasonable extrapolation of the radio emission. This point is 
well illustrated in Figure 6, in which we present six exam- 
ples of SEDs from the sources listed in Table 2. The SEDs 
of 4C41.17, BR1202-0725 and H1413+117 represent good 
examples of sources whose observed sub-millimetre SED is 
well determined and for which there can be very little doubt 
that thermal emission from dust has been detected, whilst 
the remaining SEDs illustrate the importance but also the 
limitations of a detection at a single sub-millimetre or mil- 
limetre wavelength. For example, our 800/im detection of 
53W002, if real, appears to lie clearly above the extrapo- 
lation of its radio spectrum, whereas the millimetre detec- 
tions of B2 0902+34 (Chini & Kriigel 1994) could simply be 
a detection of the high-frequency tail of the radio emission 
from the unresolved core (Downes et al. 1996) which has a 
flatter spectrum and is more luminous (with respect to the 
submillimetre emission) than, for example, the radio-core in 
the ultra-steep spectrum radio galaxy 4C41.17 (Carilli et al. 
1994). 

Figure 6 demonstrates that with appropriate target se- 
lection (e.g. ultra-steep spectrum radio sources) and with 
sensitive sub-millimetre and millimetre photometry one can 
prove that thermal emission from dust has indeed been de- 
tected, but that this must be determined with care on a 
source-by-source basis. However, even if it can be convinc- 
ingly demonstrated that the detected emission is due to 
dust, one must then address the issue of how confident one 
can be that the dust lies at the redshift of the source. We 
therefore now address this often over-looked but important 
issue. 



3.2 Potential Problem 2: What dust has been 

detected?: possible contamination by Galactic 
Cirrus 

Even if, as in the case of 4C41.17 and BR1202-0725, it 
can be shown that a rising sub-millimetre SED has been 
detected, it must be remembered that extragalactic obser- 
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Table 2. Summary of all published continuum data measuring thermal emission from dust in galaxies at z > 2 with at least one detection 
at wavelengths between 350/irra - 1.25 mm. Flux density limits are quoted at the 3<r level. 



Source name 


z 


type 


$\ ,25mm 

(mJy) 


Sg00fJ,m 

(mJy) 


S450fim 

(mJy) 


SssO/im ref. 
(mJy) 


IRAS10214+4724 


2.286 


IRAS 


24 ± 5 


50 ± 5 


273 ± 45 


1 


53W002 


2.390 


RG 




6.9 ± 2.3 




2 


H1413+117 


2.546 


BALQSO 


18.2 ± 2.0 


66 ± 7 


224 ± 38 


189 ± 56 2,3 


Q1017+1055 


3.15 


RQQ 


3.7 ± 1.2 






4 


B2 0902+34 


3.391 


RG 


3.1 ± 0.6 


< 14 


< 99 


2,5 


4C 41.17 


3.800 


RG 


2.5 ± 0.4 


17.4 ± 3.1 


< 56 


5,6 


PC2047+0123 


3.800 


RQQ 


1.9 ± 0.5 






7 


BR1117-1329 


3.96 


RQQ 


4.09 ± 0.81 






4 


BR1 144-0723 


4.14 


RQQ 


5.85 ± 1.03 






4 


8C1435+643 


4.252 


RG 


2.6 ± 0.4 


< 13 




2,7 


BRI1335-0417 


4.40 


RQQ 


10.26 ± 1.04 






4 


BRI0952-0115 


4.43 


RQQ 


2.78 ± 0.63 






4 


BR1033-0327 


4.51 


RQQ 


3.45 ± 0.65 


12 ± 4 




4,8,9 


BR1202-0725 


4.69 


RQQ 


10.5 ± 1.5 


50 ± 7 


92 ± 38 


4,8,9 



1. Rowan-Robinson et al. 1993; 2. this paper; 3. Barvainis et al. 1995; 4. Omont et al. 1996; 5. Dunlop et al. 1994; 6. Chini & Kriigel 

1994; 7. Ivison 1995; 8. Isaak et al. 1994; 9. McMahon et al. 1994 




Figure 5. A comparison of the measured SOOfim flux densities of various high-z and low-z active galaxies; starbursts (solid-stars), 
Scyfcrts (dottcd-circles), ULIRGS (crosscd-circles), radio galaxies (open-stars), radio-quiet quasars (solid-circles), taken from Hughes, 
Ward & Davies (1997), Hughes et al. (1993) and referenced data included in Table 2 (this paper). The flux densities of IRAS10214+4724 
and the Cloverleaf quasar (H1413+117) have been corrected for their amplification due to lensing, see §4. The error-bar associated with 
IRAS10214+4724 represents the uncertainty in the amplification factor at FIR wavelengths. An expanded version of this figure between 
redshifts 0.5 — » 10 is shown in the right-hand panel. The solid curve represents the predicted dependence of the 800/im flux density with 
redshift for the nearby starburst galaxy M82. Also shown are curves for a galaxy with xlO (dashed), X100 (dashed-dotted) and X1000 
(dotted) the rest-frame FIR luminosity of M82 (L FIR = 3 X 1O 1O L ). 
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Figure 6. Optical-radio SEDs of high-z radio galaxies and radio-quiet quasars showing that, except in the case of B20902+34 (sec §3.1), 
the rest- frame submillimctrc and FIR emission is in excess of an extrapolation of the total radio continuum. The SEDs of 4C41.17 and 
B20902+34 show both the core emission and the steeper large-scale lobe emission. In each case isothermal greybody spectra, assuming 
an emissivity index j3 = 1.5, are shown for 30 K (dashed), 50 K (solid) and 70 K (dottcd-dashed). The sub-millimetre and millimetre 
data are given in Table 2. The radio, infrared and optical continuum fluxes for these objects have been taken from the following papers; 
4C41.17 - Carilli et al. (1994), Chambers et al. (1990: BR1202-0725 - Isaak et al. (1994): H1413+H7 - Barvainis & Antonucci (1996), 
Barvainis et al. (1995): 53W002 - Windhorst et al. (1991): 8C1435+643 - Lacy et al. (1994): B20902+34 - Downcs et al. (1996), Yun & 
Scoville (1996), Carilli et al. (1994), Lilly (1988). 
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vations at sub-millimetre wavelengths are made in the pres- 
ence of foreground thermal emission from dust grains in the 
galactic ISM heated by the interstellar radiation field (Low 
et al. 1984), hereafter cirrus. Indeed, given its fractal nature 
(e.g. Bazell & Desert 1988), it is possible that cirrus has sig- 
nificant structure on scales less than the typical 30" — 60" 
chop-throws employed in submillimetre observations (e.g. 
Meyer 1990, Diamond et al. 1989) and hence, since this 
foreground emission may not be completely removed with 
the usual technique of position beam-switching, it may rep- 
resent a significant source of additional noise. In the limit 
this 'confusion' noise may place a fundamental constraint 
on the useful depth of searches for high-z galaxies, and it 
is therefore important to quantify the contribution of cirrus 
emission both to the current faint submillimetre observa- 
tions described in this paper, and to future submillimetre 
high-redshift surveys. 

To achieve this we have reduced the IRAS raw de- 
tector data (CRDD) at 12, 25, 60 and 100/im at the po- 
sitions of all the high-redshift objects listed in Tables 1 
& 2, using the procedure described by Hughes, Appleton 
& Schombert (1991). Not surprisingly, none of the high- 
redshift objects were actually detected. The 100/im surface 
brightness (iioo^m) of the foreground cirrus was measured 
towards all of the high-redshift galaxies after the subtrac- 
tion of any zodiacal emission which can still be significant 
in the FIR at low ecliptic latitudes. 

The average noise level at 100/im due to cirrus, on scales 
of a few arcminutes, was ~ 0.3MJy/sr which, when extrap- 
olated to 800/im, assuming an isothermal dust spectrum 
and emissivity index (3 = 2 (Rowan- Robinson 1992), gives a 
brightness of < 0.3MJy/sr, equivalent to < lmjy/beam, for 
cirrus temperatures > 12K. In other words it is extremely 
unlikely that the 800/im detections, > 7mJy, in Table 2 are 
due to cirrus unless the dust grains are significantly colder 
than the generally accepted range of cirrus temperatures, 
15K < T < 40K, (Terebey & Fich 1986, Rowan-Robinson 
1986, Van Steenberg & Schull 1988) regardless of grain-size 
or composition. However, as a caveat we mention firstly that 
if (3 7^ 2, the submillimetre emission can increase by a factor 
~ 10 as f3 varies between 2 — > 1 and, secondly, that the cir- 
rus temperatures are generally derived from a 60/im/ 100/im, 
intensity ratio which is highly sensitive to the contribution 
from a small fraction of hot grains in the cloud, and hence 
the average grain temperature that dominates the mass of 
the cloud is often over-estimated. At submillimetre wave- 
lengths we are significantly more sensitive to emission from 
a colder component than in the FIR. 

A noise level of ~ lmjy at ~ 800/im, is similar to 
the theoretical predictions of the cirrus confusion on scales 
equivalent to the chop-throw at submm wavelengths (Helou 
& Beichman 1991, Gautier, Boulanger, & Puget 1992) as- 
suming that the spatial power spectrum of cirrus, as de- 
termined at 100/im by IRAS (between scales of 8 degs. to 
2 arcmins.), continues unbroken to smaller spatial scales 
(~ 30"), and that the empirical relationship between the 
amplitude of the variations and the cirrus surface brightness 
in the FIR is valid at longer sub-millimetre wavelengths. 

While these calculations provide reassurance that de- 
tections at the level of those described in the present paper 
are unlikely to be seriously contaminated by cirrus emission, 
they also demonstrate that cirrus contamination is likely to 
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Figure 7. Surface brightness at 800/im. of galactic cirrus as 
a function of dust temperature. The curves are extrapolations 
of optically-thin, isothermal models with a grain emissivity in- 
dex /3 = 2, assuming a 100/tra intensity of lOMJy/ster (solid 
line), lMJy/ster (dashed-line), O.lMJy/ster (dot-dashed line). 
The shaded region represents the range of the surface bright- 
ness of the 800/tm detections in Table 2 and illustrates that for 
T S> 15K the contribution from cirrus is negligible. The horizon- 
tal line represents the intensity level at ~ 800/im that can be 
reached in 3 hours integration with the future generation of sub- 
millimetre receivers, and indicates that deep cosmological studies 
may suffer from source confusion due to galactic cirrus emission. 



become a more serious issue with the increased instrumen- 
tal sensitivity offered from the new generation of submil- 
limetre receivers (e.g. SCUBA). With such instruments it 
will be possible to reach noise levels of ~ 0.3mJy/beam or 
0.08 MJy/sr in just 3 hours integration and, as a result, fu- 
ture observations of high-redshift galaxies will be sensitive 
to faint, hence possibly unresolved, clumps of cirrus with 
temperatures > 15K. This result is generalised in figure 
7 which shows the surface brightness at 800/im of galac- 
tic cirrus as a function of temperature. The three extrapo- 
lated models cover the brightness range of cirrus at 100/im 
(0.1 MJy/sr < hoo^m < 10.0 MJy/sr) towards the posi- 
tions of the high-redshift galaxies given in Tables 1 & 2 and 
are typical of the cirrus background away from the galactic 
plane (b > 10°) in all but the darkest regions (e.g. Lockman 
et al. 1996). 



3.3 Potential Problem 3: Uncertainties in the 

properties and temperature of the dust grains 

If it is assumed the submillimetre continuum (A res t > 
200/im) is due to optically-thin emission from heated dust 
grains with no additional contribution from bremsstrahlung 
or synchrotron radiation, a measure of the dust mass Md 
can be determined directly from the relationship, 
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d l + zk d est B(v rest ,T) { ' 

where z is the redshift of the source, Sobs is the observed 
flux density, k r d eat is the rest-frequency mass absorption co- 
efficient, B(v rest , T) is the rest-frequency value of the Planck 
function from dust grains radiating at a temperature T, and 
Dl is the luminosity distance (see section 3.4). 

Thus, for a given cosmology and excluding the mea- 
surement uncertainty of the continuum flux density, the ro- 
bustness of dust mass determinations from submillimetre 
photometry depends on the uncertainty in kd and T. The 
existing constraints on these two parameters are now briefly 
discussed in turn. 

3.3.1 Uncertainty in kd{\) 

Draine (1990) has pointed out that the acknowledged ad- 
vantage of using optically-thin submillimetre emission to 
determine the dust mass in galaxies is offset somewhat by 
increased uncertainty in the properties of interstellar dust, 
and hence the value of k d as A is increased from FIR to sub- 
millimetre wavelengths. A reasonable estimate of the maxi- 
mum fractional uncertainty in kd at 800/im is ~ 7, with the 
values of kd(800firn) ranging between 0.04 m 2 kg -1 (Draine 
& Lee 1984) and 0.3 m 2 kg" 1 (Mathis & Whifffen 1989) with 
intermediates values of 0.15 m 2 kg -1 (Hildebrand 1983) and 
0.12m 2 % _1 (Chini, Krugel & Kreysa 1986). In order that 
we can extrapolate k d to shorter rest- wavelengths, appropri- 
ate for high-z galaxies, we have assumed that k d <x A and 
have adopted an average value of fc d (800pim) = 0.15 ± 0.09 
m 2 kg -1 . A different choice of kd should therefore only be 
expected to result in dust mass estimates which differ from 
those which we calculate in section 4, by at most a factor of 
~ 2. 

3.3.2 Uncertainty in T 

The temperature of dust grains that radiate at submillime- 
tre wavelengths and dominate the FIR luminosity in nearby 
starburst galaxies, low-metallicity dwarf galaxies, ULIRGs, 
Seyferts and radio-quiet quasars is typically 50 ±20 K (Chini 
et al. 1989a, Chini et al. 1989b, Hughes et al. 1993, Hughes, 
Davies & Ward 1997). There is no a priori reason to believe 
that the dust radiating at FIR rest wavelengths in high- 
redshift radio galaxies and radio-quiet quasars should be sig- 
nificantly different if the star formation processes that exist 
at early epochs are similar to those in the local universe. 

The fractional uncertainty in the dust mass which re- 
sults from ignorance of the dust temperature within the 
range Ti < T < T2 is given by 

Mi _ e hVr " t/kTl - 1 

Ml ~~ e h ^t/kT 2 _ I ( 3 ) 

which increases rapidly as v r est moves above the Rayleigh- 
Jeans tail of the thermal dust emission. Dust mass estimates 
derived from rest-frame sub- millimetre photometry therefore 
have the benefit of being relatively insensitive to uncertain- 
ties in temperature, as compared to dust masses derived 
from far- infrared data (i.e. Xrest < 200/im). 

At z ~ 0, where 800/im photometry samples the 
Ray leigh- Jeans tail of the thermal dust emission for Td S> 
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Figure 8. Temperature-dependent uncertainty in dust masses 
deduced from sub-millimetre photometry in various wavebands 
plotted as a function of source redshift. If we consider a fractional 
error in dust mass of greater than 5 to be unacceptable, then 
photometry at 800/im is of little use for observations of galaxies 
at z > 4, while photometry at ~ 400/tm is of little worth (on its 
own) for z > 1. 



20 K, the fractional uncertainty in the estimated dust mass 
is proportional to the fractional uncertainty in the dust 
temperature, but as figure 8 illustrates, the uncertainty in 
Mdust increases sharply with increasing redshift unless the 
observing frequency is reduced appropriately. For example, 
at z ~ 3 the difference in dust mass between assuming 
T = 30 K and T = 70 K is a factor of 15 and 5 at 450pm and 
800pm respectively. Consequently our ignorance of the most 
appropriate value of Tdust in high-redshift galaxies is a dom- 
inant source of uncertainty in the derived dust masses and it 
is worth considering the implications of such an uncertainty 
in future cosmological studies at submillimetre wavelengths. 

It is clear that the more one benefits from the increased 
flux levels offered by intercepting the thermal spectrum 
nearer the rest-frame 100pm peak, the greater uncertainty 
there is in the derived dust mass determined from a single 
submillimetre continuum measurement. As we have sum- 
marised in Table 2, almost all photometric observations at 
millimetre and submillimetre wavelengths of high-z galax- 
ies have been restricted to filters at A > 800/im, with the 
exception of IRAS10214+4724, 4C41.17, H1413+117 and 
BR1202— 0725. Thus the majority of data have been made 
in the rest-frame Rayleigh- Jeans regime (for T > 30 K) and 
consequently, despite the inability of the data to constrain 
the temperature, the uncertainty in the mass due to this 
effect alone is no worse than a factor of 5 for galaxies at a 
redshift z < 4. 

Figure 8 clearly demonstrates the value of a general phi- 
losophy when determining the dust mass from a single pho- 
tometric detection, that is to use the longest available sub- 
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millimetre or millimetre wavelength particularly for galaxies 
at z > 1. 

3.4 Potential Problem 4: Uncertain cosmology 

The luminosity distance, Dl, which appears in equation 2 is 
given by 

Dl = jB% ( Qoz + (0 ° " 2)m ° z + 1)1/2 " (4) 

The (apparently now decreasing) uncertainty in the 
value of the Hubble constant Ho obviously affects the calcu- 
lated values of dust mass to some extent, but comparison of 
the properties of low and high-redshift galaxies is unaffected 
because M dust simply scales as H^ 2 . 

More serious is our current ignorance of fio. At z ~ 4 if 
one adopts a low density universe (i.e. fio — 0.1) the derived 
dust masses are a factor ~ 5 greater than if one assumes an 
Einstein de-Sitter Universe, and this uncertainty factor rises 
to ~ 15 at redshifts z ~ 10. The value of pointing out the 
size of this uncertainty is that it puts into perspective the 
relatively minor errors which are introduced by our igno- 
rance of the precise value of such parameters as ka- 

Bearing these uncertainties in mind, we now proceed to 
calculate our best estimates of the dust masses of the high- 
redshift objects listed in Table 2, and hence to infer their 
evolutionary status. 



4 DERIVATION OF PHYSICAL 

PARAMETERS FOR HIGH-REDSHIFT 
GALAXIES 

4.1 The dust masses of high-redshift galaxies 

At first sight the preceding section may read as a rather 
depressing litany of ever-increasing sources of error. How- 
ever, our objective is not to denigrate the usefulness or im- 
portance of sub-millimetre cosmology, but rather to delin- 
eate and quantify the main sources of error, and to identify 
what strategies must be adopted to maximise the usefulness 
of sub-millimetre and millimetre continuum photometry. In 
fact, for most of the sources listed in Table 2, the current 
observational constraints indicate that we can still derive re- 
alistic and meaningful estimates of their dust masses, with 
appropriate choices of parameters. First, as illustrated in 
figure 6, for most of these sources we can be confident that 
the detected sub-millimetre emission is produced by dust. 
Second, while it appears that cirrus confusion may be a sig- 
nificant problem for deeper surveys, our estimates of the 
background cirrus noise for the sources detected here indi- 
cate that these detections are unlikely to be significantly 
effected. Third, we can minimize the effect of disagreement 
of kd by choosing an average value. Fourth, for at least some 
sources (e.g. 4C41.17, BR1202-0725, H1413+117) there ex- 
ists sufficient observational data to set a meaningful upper 
limit on the dust temperature, and hence lower-limit on the 
dust mass. Fifth, while it is obviously beyond the scope of 
this paper to address the issue of the true value of fio, by 
adopting a high-density universe we can at least ensure that 
our dust mass estimates are conservative. 

We have therefore calculated dust masses for the sources 



listed in Table 2 using equation 2, assuming kd(800fim) = 
0.15m 2 kg- 1 , (3 = 1.5, T dust = 50K and fi = 1. The results 
are presented in Table 3 and are in the range 1 x 10 s — 2 x 
10 9 Mq allowing for amplifications of order ~ 11 (Barvai- 
nis et al. 1994) and ~ 10 - 30 (Eisenhardt et al. 1996) for 
the FIR luminosity in the lensed sources H1413+117 and 
IRAS10214+4724 respectively. 

These dust mass estimates can be used to infer the evo- 
lutionary status of the galaxies in two alternative ways. First 
one can estimate the 'current' star- formation rate of the 
galaxy. Second one can attempt to estimate the amount of 
molecular gas in the galaxy which remains to be converted 
into stars at the epoch of observation. We now consider each 
of these calculations in turn. 

4.2 Star formation rates in high-redshift galaxies 

As described in the introduction, one can estimate the 
star-formation rate at the time of observation from their 
rest-frame FIR luminosities which have been calculated be- 
tween rest-wavelengths of 2 mm — ► lO^im assuming an 
optically-thin isothermal 50K greybody spectrum (see Ta- 
bles), which can be reasonably justified in the case of 
4C41.17, BR1202-0725 and H1413+117 (figure 6, see also 
§4.2.1), normalised to the millimetre and submillimetre flux 
densities in Table 2. In the absence of any contribution from 
an AGN or amplification due to lensing, the rest-frame FIR 
luminosities, which lie in the range 4 x 10 12 — > 6 x 10 13 Lq, 
imply large current SFR's of greater than several 100M Q /yr. 
A comparison of the 1000 — 8^im luminosities of low-z 
ULIRGs (Sanders et al. 1991) with those calculated from 
a single 50 K isothermal model indicates that we have un- 
derestimated the rest-frame FIR luminosities of the high-z 
galaxies by a factor of ~ 2.5 if they contain contributions 
from hotter dust (T > 100K) components which radiate 
strongly at mid-IR wavelengths. Hence the FIR luminosities 
and SFRs in Table 3 should be treated as lower-limits when 
compared to those of lower-redshift galaxies. Further, albeit 
controversial evidence for extreme SFRs and young galaxy 
ages (< 1 Gyr) have been found in the rest-frame UV-optical 
morphologies and SEDs of 53W002 and 4C41.17 (Windhorst 
et al. 1991, Chambers et al. 1990, Mazzei & de Zotti 1995). 

Thus, despite the uncertainties described in §3 we can 
still conclude that the high-z radio galaxies and radio-quiet 
quasars, which have been detected at submm-mm wave- 
lengths, are extremely dusty, with dust masses > 10 x larger 
than observed in their low-z (z < 0.5) counterparts (Chini 
et al. 1989a, Knapp & Patten 1991, Hughes et al. 1993), and 
with levels of starformation and starforming efficiencies sim- 
ilar to, or exceeding those observed in low-z ULIRGS (see 
figure 9) . 

Particularly if regarded as lower limits, these inferred 
SFRs arc sufficiently high to be at least consistent with the 
formation of a giant elliptical galaxy in ~ 1 Gyr. However, 
such a calculation is open to two important criticisms. First, 
for the active objects considered here it is unclear whether 
the dust is heated primarily by young stars or by a quasar 
nucleus, although this argument can be at least partially 
countered by the fact that the very existence of such large 
masses of dust indicates extensive recent star-formation. A 
second, and more important criticism is that an estimate of 
the 'current' SFR can never shed any light on the length of 
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Table 3. Dust masses, molecular gas masses, FIR luminosities and starformation rates (SFRs) of high-z galaxies. Corrections for the 
amplification of the FIR emission due to lensing have been applied to the values for IRAS10214+4724 and H1413+117 (see §4.1). 
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time over which such a high star-formation rate has been 
maintained. In other words it is clearly impossible to dis- 
tinguish a brief, albeit violent burst of star-formation ac- 
tivity (triggered, perhaps, by an interaction) from the sort 
of sustained high-level star-formation required to produce 
~ 10 11 1O 12 M of stars. What is required is some es- 
timate of the past and future star-formation history of the 
galaxy. Such an estimate can be provided, to first order, by 
using the calculated dust mass of a galaxy to estimate its 
gas mass. 



4-2.1 Are simple isothermal models adequate? 

We have assumed that the rest-frame FIR luminosities in 
table 3 are dominated by radiation from grains with an av- 
erage temperature of 50K, and that the SFRs can be derived 
from a single isothermal component. This is of course a sim- 
plification, since real galaxies clearly have a distribution of 
grain temperatures that contribute to the rest-frame emis- 
sion between 200/im - 50/im. However the thermal emission 
from grains which radiate at a temperature T, with an ef- 
ficiency Qabs oc v 13 , has a peak in the spectrum of Si, at a 
wavelength \(pm) ~ t(k) (3T/3)' anc ^ consequently there 
remains only a restricted range of temperatures (20-70 K) 
that contribute significantly to the the rest-frame FIR. 

In figure 10 we quantify the discrepancy between our 
calculations of the FIR luminosities in table 3 and, where 
overlap exists, the bolometric starburst luminosities calcu- 
lated from the more sophisticated radiative transfer models 
of Rowan- Robinson (1996) and the evolutionary synthesis 
models of Mazzei & de Zotti (1996). Figure 10 demonstrates 
that in the absence of sufficient observational data to con- 
strain models between 800/i — > 2p, the adoption of a single 
average grain temperature (50K) provides extremely good 
agreement, to within 20%, between the FIR and bolometric 



luminosities (and also the current SFRs) for galaxies ex- 
pected to be undergoing vigorous starformation. 

4.3 Molecular gas masses in high-redshift galaxies 

A measure of the total (HI + H2) gas mass of a high-redshift 
galaxy is a very useful indicator of its evolutionary status 
because, as long as one can be reasonably confident that the 
galaxy in question is the progenitor of a present-day giant 
elliptical (as is the case for radio galaxies) , it provides a mea- 
sure of the fraction of the galaxy which has yet to be turned 
into stars at the epoch of observation. One way to determine 
the gas mass of a galaxy is to estimate its molecular (-H2) 
mass from a measurement of the intensity of CO. Unfortu- 
nately, in the absence of lensing, it has proved impossible to 
detect the molecular gas in high-z galaxies directly (van Ojik 
et al. 1997, Evans et al. 1996, Barvainis & Antonucci 1996), 
with the possible exception of BR1202— 0725 (Omont et al. 
1996, Ohta et al. 1996). However, an alternative and poten- 
tially more productive approach is to convert dust mass to 
gas mass by adopting a 'reasonable' value for the gasidust 
ratio, Mif 2 /M<j. However Mh 2 /M<j is not a well determined 
quantity in galaxies in the local universe, let alone at high 
redshift. Studies of damped Lyman-a systems (DLAAS) cur- 
rently provide the only opportunity to directly measure the 
dust content of the universe at early epochs. At z ~ 3 it has 
been suggested that the gas-dust ratio in DLAAS is 400- 
2000 (Fall, Pei & McMahon 1989, Pettini et al. 1994), a 
value significantly higher than the galactic value of 100-160 
(Hildebrand 1983, Savage & Mathis 1979), and on average 
higher than that found in nearby spirals (~ 500, Devereux 
& Young 1990), ellipticals (~ 700, Wilkind et al. 1995) and 
ULIRGS (540 ± 290, Sanders et al. 1991). 

While it may be true that the ratio Mn^jMd in high- 
redshift Lyman-a absorbers is significantly greater than in 
present day galaxies, it seems likely that high-redshift ra- 
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Figure 9. A representation of the data given in Table 3 for the detected high-z radio galaxies and radio-quiet quasars, compared to 
their lower-redshift counterparts, ULIRGs and elliptical galaxies. The symbols have the same meaning as in figure 5, with additional 
low-z powerful radio galaxies (Mazzarella et al. 1993), ULIRGs (Sanders et al. 1991) and elliptical galaxies (Lees et al. 1991) detected 
in CO marked as solid stars, crossed-circles and solid squares respectively. The diagonal lines indicate constant Lpir/WIh 2 , while the 
vertical dashed-line shows the gas mass boundary, to the right of which one can expect to find the projenitors of the most massive 
elliptical galaxies. The parallelogram encloses the high-z galaxies, which lie along the lower boundary, and depicts the typical range of 
the increase in their FIR luminosities and SFRs due to a contribution from hotter (T > 50K) dust (see §4.2). This figure demonstrates 
that it is difficult, except possibly in the case of BR1202— 0725 (although it may be lensed, see §5), to describe any of the high-z galaxies 
detected at submm and mm wavelengths as genuinely primeEval, since they lie outside the shaded region marked " PGs" which represents 
the parameter space populated by primaeval galaxies, with our definition of Mj/ 2 > 5 X 10 11 Mq and a SFR> 500 Mq. Varying the 
assumed dust temperature through a reasonable range (30-K" — » 70 K), as illustrated by the representative locus passing through 4C41.17, 
struggles to change this basic conclusion, except perhaps for 4C41.17 itself. 



dio galaxies, being the progenitors of present day ellipticals, 
are considerably more evolved than Lyman-a absorbers at 
comparable redshift (which, it has been argued, are the pro- 
genitors of present-day disc galaxies). Thus there seems to 
be no clear justification for adopting an extreme gas:dust 
ratio and so we have decided to adopt a conservative value 
of MH 2 /Md ~ 500, consistent with the values reported for 
galaxies in the local universe (the appropriateness of this 
chosen value is discussed further in the next section). 

The resulting values of Mh 2 determined from contin- 
uum measurements are listed in Table 3, and, are < 8 x 
10 11 Mq. These estimates of the H2 gas mass limits are con- 
sistent with the failure to detect large reservoirs of molecular 
gas directly in various CO transitions (Evans et al. 1996; van 
Ojik et al. 1997, Barvainis & Antonucci 1996) a fact which 
lends some circumstantial support to our adopted gas:dust 
ratio. Nevertheless, these gas masses are still extremely large 
and represent a significant fraction of the present day stel- 



lar mass of the largest elliptical galaxies. The crucial issue of 
whether this 'significant fraction' could be sufficiently large 
to indicate that these objects deserve to be described as 
'primaeval' is considered below in the final section of this 
paper. 



5 CONCLUSIONS: THE EVOLUTIONARY 
STATUS OF HIGH-REDSHIFT GALAXIES 

We have presented sensitive continuum measurements at 
800/im of a sample of high-redshift radio galaxies and 
quasars. These observations were motivated by the goal of 
determining their evolutionary status in a relatively model 
independent manner via a measure of their dust masses, 
star-formation rates, and molecular gas masses. In addition 
to our 800/xm data presented here, we have taken the op- 
portunity to gather together all other published millimetre 
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Figure 10. A comparison of the rest-frame FIR luminosities 
(Lpjn) calculated from a 50K isothermal model (table 3) and the 
bolomctric luminosities (L(, ;) determined from radiative trans- 
fer models of Rowan-Robinson (1996, MRR) - filled squares, and 
from the range of evolutionary synthesis models of Mazzei & de 
Zotti (1996, MdZ) - lines, which provide a reasonable fit to the 
submm, optical and IR data. The shaded region represents a max- 
imum discrepancy of 20% about a perfect correlation between the 
FIR luminosities, and hence SFRs, calculated from isothermal 
and more sophisticated models. 

and submillimetre observations of high-redshift galaxies and 
quasars at z > 2, and Table 3, in effect, summarizes the cur- 
rent status of sub-millimetre cosmology, 

While it is expected that this subject will be revolution- 
ized by the significantly more sensitive observations which 
should be possible in the very near future using new submil- 
limetre continuum bolometer arrays such as SCUBA (Gear 
& Cunningham 1995), we have argued in this paper that 
there is much to be learned from considering the uncertain- 
ties which afflict the reliable interpretation of the existing, 
albeit sparse dataset. In this spirit we conclude by consid- 
ering what can be deduced about the evolutionary status 
of high-redshift galaxies from the numbers presented in Ta- 
ble 3. 

To focus the argument we adopt the strict (and de- 
liberately extreme) definition that a genuinely young el- 
liptical galaxy would be expected to display a very high 
star-formation rate (> 500M Q j/r _1 ), and to contain a to- 
tal HI + H2 gas mass equivalent to the majority of the 
stellar mass found in the most massive present-day giant el- 
lipticals (e.g. Mhi+h 2 > 1O 12 M0). The adoption of such a 
large present-day mass can be justified on two levels. First, if 
one wants to unambiguously prove that a given high-redshift 
object is primaeval on the basis of its gas mass, one needs 
to demonstrate that, even if its eventual present-day stellar 
mass is extremely large, the bulk of it remains in gaseous 
form at the high-2 epoch of observation. Second, since it 
is now well-established that the most powerful radio sources 



and quasars in the low-redshift universe reside in host galax- 
ies with luminosities several times greater than L* (Taylor 
et al. 1996), it is not unreasonable to assume that the con- 
siderably more luminous active objects listed in Table 3 lie 
in still deeper potential wells. 

Since we are unable to detect the contribution from 
21 cm HI emission, which is shifted to 250-500MHz at high- 
redshift, we use the mean ratio of Mhi/Mh 2 ~ 1-0 ± 0.9 
found in elliptical, early-type and interacting galaxies (Lees 
et al. 1991, Young & Knezek 1989), to place a lower mass 
limit of Mh 2 > 5 x 1O 11 M0 on the molecular gas content of 
the most luminous primaeval or proto-elliptical galaxies. 

As illustrated in Figure 9, when judged against these cri- 
teria of SFR and Mh 2 , although several of the detected ob- 
jects in Table 3 may perhaps have adequate star-formation 
rates, only one source, the quasar BR1202— 0725, lies within 
the parameter space populated by primaeval galaxies under 
the above definition. Given that both IRAS10214+4724 and 
the Cloverleaf quasar would also have qualified as genuinely 
young on the basis of the H2 gas mass prior to correcting 
for the estimated effects of lensing, the suspicion, borne out 
by the recent CO detection (Omont et al. 1996, Ohta et 
al. 1996) must be that BR1202— 0725 may also have am- 
plified emission. Therefore based on our analysis, the sub- 
millimetre properties of high-redshift AGN, which have been 
identified at optical/IR/radio wavelengths, are more typi- 
cal of a strong starburst in an otherwise well-formed galaxy 
than of genuinely primaeval objects. Such a result could in- 
dicate that we are seeing the final stages of formation of 
these massive objects, or a starburst triggered by a galaxy- 
galaxy interaction which is also responsible for the observed 
AGN activity (in a heirarchical picture of galaxy formation 
these two alternatives could essentially amount to the same 
thing). 

In judging the robustness of the basic conclusion that 
these high-z galaxies are not primaeval, it is important to 
consider in turn the effect of the dust temperature and 
adopted gas:dust ratio on the dust masses (— ► Mh 2 ) and 
FIR luminosities (— > SFR) in Table 3 since both of the lat- 
ter quantities are determined from an isothermal grey-body 
fit to the submillimetre continuum. 

First, the combined effect of the temperature uncer- 
tainty (T d ~ 50 ± 20K) on the FIR luminosities and dust 
masses, which was discussed at length in §3.3.2, is illus- 
trated in the representative locus passing through 4C41.17 
in figure 9. This locus, which is appropriate for galaxies at 
z — 2 — 5 shows that the upper-limit in temperature (70K) 
moves the high-z sources away from the parameter space oc- 
cupied by primaeval galaxies and into a region requiring ex- 
tremely high star-forming efficiencies (S> 100Lq/Mq). Con- 
versely, while reducing the temperature to 30K can increase 
the inferred gas mass to close to the required value, this is at 
the expense of the FIR luminosity and star-formation rate 
which is reduced to < 500MQjyr _1 . In summary, our basic 
conclusion - that we are not seeing the primary formation 
event of these galaxies - is, with the possible exception of 
4C41.17, unaffected by varying the temperature of the dust 
between reasonable limits. 

It is worth noting that while our adoption of a dust 
temperature of 50 K for these high-redshift objects results in 
values for Lfir/Mh 2 — 100 somewhat greater than is seen 
in isolated systems at low-redshift (12 ± ZLq/Mq) (Young 
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et al. 1986, Solomon & Sage 1988), such values are not un- 
reasonable and are comparable to that seen in the most lu- 
minous low-redshift interacting galaxies (78 ± 14L©/M©). 
Indeed the location of the high-redshift sources on Figure 9 
seems perfectly consistent with an extrapolation of the up- 
per envelope of star-forming efficiency defined by the most 
luminous ULIRGs. 

Second, we address the issue of whether the uncer- 
tainty in the assumed ratio Mn 2 /Md ~ 500 (§4.3) can affect 
our conclusions. Obviously it can if a 'reasonable' range of 
Mn 2 /M d at high-redshift is taken to extend to values of 
several thousand as suggested by the analysis of Fall et al. 
(1989). However, as already discussed above, the relevance of 
the investigations of Lyman a absorbers to the present study 
of high-redshift AGN is dubious. More significant is the 
good agreement between our adopted value and the mean 
MH 2 /M d ratio, (426 ± 96), measured in the only 3 high- 
z sources (H1413+117, IRAS 102 14+4724, BR1202-0725) 
with confirmed CO line and submillimetre continuum de- 
tections. The gasidust ratio in these high-z AGN should be 
relatively unaffected by lensing assuming that, as in low-z 
galaxies, the CO and submillimetre continuum have simi- 
lar spatial distributions, and hence similar amplifications. 
Thus, taken at face value, these few CO detections of high- 
redshift objects support our view that MH 2 /Md ~ 500 is 
appropriate for luminous galaxies out to z ~ 5, and ex- 
plains why sub-millimetre continuum observations of high- 
redshift objects have been more successful than attempts 
to detect the molecular gas directly through CO line ob- 
servations (and casts even more doubt on the claimed CO 
detection of 53W002, Yamada et al. 1995). 

To summarize, figure 9 demonstrates that adoption of 
a dust temperature of 50 K and a gas:dust ratio of 500 leads 
one to conclude that the high-redshift objects listed in Ta- 
ble 9 would be better described as highly efficient starburst 
galaxies than genuinely primaeval objects. This conclusion is 
relatively immune to alteration of the adopted dust temper- 
ature. Indeed, to alter it significantly requires the adoption 
of a gas:dust ratio ~ 5 times greater than has been assumed, 
and the low success rate of molecular CO line observations 
argues against this option. However, to finish on a caution- 
ary note, there is one way to increase the inferred gas mass 
without increasing its detectability through line observa- 
tions, and that is to relax the assumption of an Einstein-de- 
Sitter Universe. Ignorance of Qo is of course a problem which 
afflicts studies of galaxy evolution at all wavelengths but, as 
discussed in section 3.4, at z ~ 4 adoption of Qo ^0.1 
instead of Qq = 1 increases the inferred dust mass by a fac- 
tor of 5 and so, all other things being equal, would increase 
the values of both Lfir and Mh 2 plotted in Figure 9 by 
a factor of 5 for the high-redshift objects only. Under these 



circumstances, with star-formation rates (~ 10 M s yr" ) 
and molecular gas masses (Mh 2 — 5 x 10 11 M Q ) it would 
be hard to escape the conclusion that essentially all of the 
high-redshift objects so far detected at sub- millimetre wave- 
lengths are in fact primaeval galaxies (one could, for exam- 
ple, reduce the inferred gas mass by adopting a high dust 
temperature such as T ~ 70K, but this would in turn imply 
a truly prodigious star-formation rate). 

Thus the unambiguous determination of the evolution- 
ary status of the high-redshift objects detected to date at 
sub- millimetre wavelengths must await improved constraints 



on fio- However, we regard it as encouraging for the future 
of sub-millimetre cosmology that the uncertainties peculiar 
to the interpretation of these sub-millimetre data can be 
justified as being comparable in size to the current cosmo- 
logical uncertainties. In general, continuum observations at 
sub-millimetre and millimetre wavelengths have proven to 
be more successful and have provided a comparable or lower 
H2 mass limit than spectral line data obtained in an equal 
integration time. This observing efficiency, together with 
an expected improvement in our understanding of the dust 
emissivity and gas-to-dust ratio at high-redshift, justifies the 
continuation of this approach with the future generation of 
sub- millimetre continuum instruments (e.g. SCUBA) which 



will be sensitive to SFRs < 100M Q yr _1 at z 



5 and 



hence will detect unlensed versions of IRAS10214+4724 and 
H1413+117 with ease and also detect the equivalent of low- 
redshift ULIRGS (e.g. Arp220, Mrk231) at z > 4. Given 
that the sub-millimetre properties of the high-redshift radio 
galaxies and quasars studied to date point towards a more 
dramatic formation event at yet higher redshift, it will be 
of particular interest to discover whether any of the objects 
detected by sub- millimetre surveys meet the criteria defined 
in Figure 9 for a genuinely primaeval giant elliptical galaxy. 
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